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Figure Captions 



Figure 1 


Naval Postgraduate School R/V Acania at the conclusion of 
CKWCOM-76 . 


Figure 2 


Normalized heat conduction coefficient (eddy diffusivity) 
versus Richardson number. 


Figure 3 


Variation of atmospheric quantities at 2 = 10 meters: 

a) temperature structure function, C T 

b) rate of dissipation of velocity fluctuations, £ 

c) wind velocity, U 


Figure 4 


2 

Temperature structure function, , at Z = 10 meters 

versus time in Monterey Bay. Solid circles are unstable 
conditions (12/76) and solid triangles are near neutral 
conditions (7/76) . 


Figure 5 


2 

Diurnal variation of C ^ : 

a) open ocean 

b) coastal ocean 


Figure 6 


R/V Acania tracks off Southern California during fog events 
for CEWCOM-76 : 

a) Sep 27-29 (Fog #1) 

b) Oct 5-9 (Fog #'s 2,3,4) 

The six digit number indicates the day of the month (dd) and 
time of day (tttt) in the code - ddtttt. 


Figure 7 


Sea surface temperature contours from aircraft IR measurements 
Courtesy of Backes (1977) and Ralph Markson, Airborne Research 
Associates . 
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Figure 8 C T ~ at 2 = 10 meters versus time for fog events: 

a) Fog #1 Sep 27-28 

b) Fog #2 Oct 5 

c) Fog #3 Oct 8 

d) Fog #4 Oct 9 



Figure 9 • £ at Z = 10 meters versus time for fog events : 

a) Fog #1 Sep 27-28 

b) Fog #2 Oct 5 

c) Fog #3 Oct 8 

d) Fog #4 Oct 9 



Figure 10 Sensible heat flux, W, versus time for fog events: 

a) Fog #1 and Fog #2 

b) Fog #3 and Fog #4 

Figure 11 Radiosonde measurements of temperature (solid line) 

and relative humidity (dashed line) versus altitude 
for fog events : 

a) Fog #1 

b) Fog #2 

c) Fog #3 

d) Fog #4 

Courtesy of NWSMET. 



Figure 12 Ensemble average of turbulence parameters one hour before 

and after fog is encountered and one hour before and after 
fog dissipates: 

a) C t 2 /C t 2 versus time 

b) e/£ versus time 
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Figure 13 



Height dependence of normalized to the Z = 6.6 meter 

value for fog (x's) and non-fog (o's) data. Neutral 
equilibrium yields the line indicated, non-neutral 

equilibrium would be a less steep line in this representa- 
tion . 

Figure 14 Dimensionless temperature structure parameter versus 

Richardson number (atmospheric stability) for fog (x's) 
and non-fog (o's) data. The solid line is the overland 
results from Wyngaard, et. al . (1971) . 
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Atmospheric Turbulence Measurements in Marine Fog 
during CEWCOM-76 

C.W. Fairall, K.L. Davidson, T.M. Houlihan and G.E. Schacher 

Environmental Physics Group 
Naval Postgraduate School 
Monterey, California 93940 

Abstract 

Atmospheric turbulence parameters £ (turbulent energy dissipation 
2 

rate) , and (temperature structure function) were measured in the 

marine boundary layer during the Cooperative Experiment for West Coast 

Oceanography and Meteorology (CEWCOM-76) . The atmospheric stability 

was characterized by measurements of profiles of mean wind velocity, 

temperature and relative humidity. This paper will focus on turbulence 

properties and turbulent heat flux measured during four marine fog 

events . The fogs were characterized by a slightly unstable boundary 

layer with moderate increases of heat flux during the fog. An ensemble 

2 

average of the fog events showed that C^ and £ reached peak values 

immediately before and after the fog was encountered. An analysis of the 

2 

dimensionless temperature structure function (C^ normalized by the square 
of the temperature gradient) gave lower values of temperature fluctuations 
in fog than in clear air. 
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I. Introduction 



During September and October of 1976 the Naval Postgraduate School 
(NPS) shipboard turbulence and profile systems were used in a marine fog 
and boundary layer research cruise off Southern California aboard the NPS 
research vessel R/V ACANIA (Fig. 1). Termed the Cooperative Experiment in 
West Coast Oceanography and Meteorology (CEWCOM-76), the cruise was jointly 
organized by NPS and the Naval Electronics Laboratory Center (NELC) of San 
Diego in collaboration with CALSPAN Corporation under the sponsorship of 
NAVAIR. Also represented on board were the Naval Avionics Facility, 
Indianapolis (NAFI) , the Naval Research Laboratory (NRL) , the California 
Air Resources Board (CARB) , and a Naval Weather Systems Service Mobile 
Environmental Team (NWSMET) . This was the most recent of a series of data 
gathering cruises forming a major part of a program of marine fog investi- 
gation sponsored by Naval Air Systems Command. It was also the third such 
cruise to include measurements of atmospheric turbulence. The first was 
conducted in 1974 aboard the R/V ACANIA off Northern California (Mack et al, 
1975); the second was conducted in 1975 aboard the USNS HAYES off Novia 
Scotia (Gathman and Larson, 1976). 

Marine fog is one of the least understood and most complicated atmos- 
pheric phenomena. Mack (1975) has identified at least five different fog 
formation "scenarios' 1 over the ocean. Although considerable fog forecasting 
success has been obtained by skilled phenomenologis ts employing various 
fog indices (such as inversion height and air sea temperature difference) 
in local coastal areas (Leipper, 1948, Backes , 1977), accurate marine fog 
prediction in the open ocean will surely benefit from applications of know- 
tecge from the synoptic to the microscale. On the microscale, the pro- 
duction of fog is governed by several (interacting) factors: aerosols, 



9 



radiation, hydrostatic stability, temperature inversion characteris tics , 
and turbulent transport of heat and water vapor. Turbulent transport was 
first incorporated into fog modeling by Rodhe (1962). More modern fog 
models (Lee and Lin, 1975; Barker, 1975) employ sophisticated turbulent 
parameter stability and scaling laws (see section II) that are unverified 
in the marine boundary layer (Davidson et al, 1977). The main thrust of 
the NPS turbulence group during CEWCOM-76 was to attack these problems by: 

1) providing modelers with measurements of turbulence and mean profiles 
in fog, 2) obtaining data for evaluation of the turbulence scaling laws 
in the marine boundary layer, 3) examining some of the characteristics of 
turbulence "background" levels for non-fog conditions in the marine boundary 
layer. 
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II. Theory 

A. Boundary layer equations 

The boundary layer is that part of the atmosphere where friction 
with and heating by the surface play an important role in the generation 
of turbulence. Near the surface the shear stress and scalar fluxes are 
essentially constant with height. In this region the fluxes can be rep- 
resented in terms of scaling parameters (such as u* and T*) that are in- 
dependent of height. We shall refer to this layer of nearly constant shear 
stress and flux as the surface layer. For a complete treatment of the 
surface layer equations we suggest Lumley and Panofsky (1964), Businger 

(1971) and Kraus (1972). The normalized fluxes of momentum (F ), heat 

m 

(F, ) and water vapor (F ) are 
h q 



F m = “ <u'v’> 


2 

= u * 


(la) 


F h = - < T'w’> 


" u * T * 


(lb) 


F = - < q'w'> 


* 

t? 

* 

d 

n 


(lc) 



where u* and w' are the fluctuating horizontal and vertical wind velocities, 

respectively, T* is the fluctuating potential temperature and q* is the 

fluctuating specific humidity. The Reynolds shear stress, T, is related to 

2 

the friction velocity, u*, by T = p u* , where P is the density of air. In 
Che surface layer we can write simplified turbulent energy balance equations 
using production terms (proportional to the vertical gradient, 3X/3Z j and 
dissipation terms [e ] as follows: 



velocity: u* 3u/3z - | = £ 

temperature: u^ T* 3T/3z = £^ 
humidity: u^ q* 3q/3z = £^ 



(2a) 

(2b) 

(2c) 
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Note that the velocity equation (where = e) has an extra pro- 
duction (or suppression) term due to the bouyancy of air (T^ is the vir- 
tual potential temperature). In the surface layer, the appropriate tur- 
bulent transport scale parameter is the height above the surface, z. In 
this regime, turbulence properties can be scaled through the similarity 
parameter z/L, where L is the Monin-Obukhov stability length. 



T u. 



L = 



<g T 



(3) 



where < = .35 is the von Karmon constant 



file functions <J> x (z/L) as (Businger,et 



. r / j \ <z dx 

<Z/L) ■ ^ 37 



We can now define the mean pro- 
al, 1971) 



(4) 



If we define a dimensionless dissipation function (Wyngaard and Cote, 1971) 



E x (z/L) 



£ KZ 
X 




(5) 



then the energy equations become 



*u 


(z/L) - z/L = E u (z/L) 


(6a) 




(z/L) = E t (z/L) 


(6b) 


<P 


(z/L) = E (z/L) 


(6c) 


q 


q 



Present estimates of 4 and E are given in Appendix A. 

X u 

In the inertial subrange, turbulence is nearly isotropic and we 
can represent the intensity of turbulence fluctuations as a one-dimensional 
power spectrum, S x (k), of the form proposed by Kolmogorov 

S x (k) = .25 C x 2 k~ 5/3 (7) 
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where k is the wavenumber and C “ is the struciure function. If the 

x 

quantity x has value x(r) at position r and value x(r + d) at position 
r + d then 



C x 2 = <^[ x(r) - x (r + d)} 2 } d 



( 8 ) 



2 

In the inertial subrange is independent of d. Corrsin (1951) has shown 

that the structure function can be related to the dissipation by 



C 



2 

x 



4 g e 1/3 e 
x x 



(9) 



where 6 is a constant, B- = B = .81 and 8 = .52. We can construct a 

x T q u 

dimensionless temperature structure parameter, f(z/L), by combining eqs. 
6 and 9 (Wyngaard, et al. , 1971) 



2 ? -?/3 

c/ = T*- z f 1 (z/L) 



( 10 ) 



or, using eq. 4 



C T 2 = (3T/3z) 2 Z 4/3 f 2 (z/L) 



(ID 



We can also use eq. 5 to obtain a similar expression for £. 
,1 3 



z = 



E (z/L) 

Kz u 



( 12 ) 



Since the Monin-Obukhov length is difficult to measure, we normally use 
an alternate stability variable-the gradient Richardson number, R , 



(31 /3z) 

R = -§ 4 o 

i T (3u/3z) *■ 



(13) 
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which is related to z/L by 



R ± = (z/L) <t> T (z/L) / ($ u (z/L)) 2 (14) 

2 

The height dependence of and £ can be obtained from eqs. 10 and 12; for 

2 -2/3 -1 

near neutral conditions (R^ - 0) - z and £ - z . For unstable 

2 -4/3 

conditions (R^ << 0) - z and £ - constant. For stable conditions 

(Rf > 0), turbulent eddies are increasingly damped by increasing stable 
stratification of the atmosphere. Overland measurements indicate complete 
damping of turbulence for R_^ > .21. 

B. Flux calculations 

There are four methods of calculating scalar fluxes: eddy correlation, 

gradient, turbulence, and bulk. The eddy correlation method is a direct 
measurement of x*w f in equation 1. This requires a stable platform so 
we cannot use it for shipboard measurements. The bulk method (Friehe, 1977) 
relates the flux to the mean wind velocity and the air-sea temperature or 
humidity difference. We have found this method to be unsatisfactory due 
to difficulties in measuring the sea surface temperature. The gradient 
method relates the vertical gradient of the scalar to the flux through the 
eddy diffusivity, For instance, the normalized heat flux is 

F h = - Kf OT/3z) (15) 

Unfortunately, K^, is a function of stability, 

*4 - < z u*/4> t (z/L) (16) 

shown in Fig. 2. Note that K^, = 0 for R^ > .21 in this model. The tur- 

2 

bulence method is based on calculations of T* from and £. Combining 

eqs. 9 and 6b we find 
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(17) 



T 



2 

* 



2 ] / 3 

< z C T £ ' 7(3.2 



u* 4 > t (z/L) ) 



Substicucing for u* with eq. 12, we obtain 



F h = 



+ r £ 
- £ T e 



1/3 



J&sl 



4/3 



1/2 



3.2 (E u (z/L)) 



1/3 



1> T (z/L) 



(18) 



The sensible heat flux is related to the normalized flux by 

W = C F h (19) 

3 3 

where C = 1.3x10 joule/ (°C m ) is the heat capacity of air. Equations 15 
and 18 have their analogous forms for calculation of water vapor flux. Un- 
fortunately, instrumental limitations and problems (see Section III) have 

2 

prevented us from measuring C and Oq/8z) during fog events. 
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III. Instrumentation and Analysis 

Since the NPS work represented only a part of the total shipboard 
effort for CEWCOM-76, we will first list the key measurements made by the 
other groups: 

1. CALSPAN - aerosol spectra and chemical composition, visibility 
and profiles. 

2 . NR L (R. Jeck) - Knollenberg aerosol spectrometer. 

3. NAFI (J. Russell) - microwave ref ractometer . 

4. CARB - ozone, carbon monoxide and NO . 

x 

5. NWSMET - twice daily radiosondes. 

The NPS instrumentation was primarily devoted to measurement of mean 
and turbulence profiles. A comprehensive description of this system is 
the subject of another paper (Houlihan, et al., 1977). The profiles were 
measured at four fully instrumented levels plus a sea surface temperature 
sensor. The level heights were 4.2, 6.6, 11.3 and 17.7 meters above the 
mean waterline. 

Mean temperature was measured with Hewlett-Packard Model HP2801A quartz 
oscillator thermometers, accurate to .01 °C. Mean humidity was measured 
with Hydrodynamic , Inc . Digital II type 15-1818 LiCl Dunmore sensors, ac- 
curate to about 3% RH. These humidity sensors have very slow response 
above about 95% RH and response variations from sensor to sensor make them 
unreliable under fog conditions. For this reason, we do not have fog 
event humidity profiles of sufficient accuracy to calculate the humidity 
flux. The wind velocity was measured with Thornthwaite Model 101 cup 
anemometers, accurate to about 5%. 

Wind speed fluctuations (for e) were measured with Thermo-Sys terns , Inc. 
Model TSI 1054B constant temperature anemometers using TSI 1210 probes with 
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6 mil. cylindrical hoc-film sensors. The hot films were calibrated on a 

TSI il25 Calibrator with Che cables actually used during the shipboard 

2 

measurements. Temperature fluctuations (for C ) were measured with GTE 

Sylvania ac wheatstone bridges (3.0 kHz carrier) using TSI 1210-P.8 micro 

thermal probes w^th 2.5 li diameter platinum wire. We had planned to measure 

2 

water vapor fluctuations (for ) at one level using an Electro-Magnetics 
Research Corporation Lyman-alpha sensor, however this instrument failed 
early in the cruise. The mean data was continuously averaged and periodi- 
cally logged by a microprocessor designated MIDAS (Microprogrammable 
Integrated Data Acquisition System) developed at NPS and adapted for the 
ACANIA by Plunkett (1977). The fluctuating data was amplified, filtered 
and recorded on a Honeywell Model 5600 FM tape recorder. 

In addition, the NPS acoustic radar (Aeroenvironment Model 300) was 
mounted on the ship and provided a continuous monitor of the inversion 
height. 

The analysis of the turbulence signals was auided by noise consider- 
ations. The combination of long sensor cables (about 150 ft) and somewhat 
impure ship generated power restricted our analysis options — particularly 
for the hot film systems. The velocity signals were processed with Fourier 
power spectrum analyzers to yield S^(f). Using Taylor’s hypothesis 
(k = 2*17 f / u ) and equations 7 and 9 we find 



e ■ 2.67 ^ [f 5/3 S u <f>] 



3/2 



( 20 ) 



The spectral method has the advantage that power source noise 

2 

(60,120,180 Hz, etc.) can be ignored. The bulk of the data were ob- 

tained from paired sensors (separated by 0.3 meters) using eq. 8. In some 

2 

cases only one sensor was operable; these were spectrum analyzed and 
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was calculated using eq. 7. Being based upon the high frequency turbulence 
components, both of these techniques are relatively unaffected by ship 
motions . 



18 



IV. Results 



A. General Background 

Before plunging into the fog data, it is of interest to examine 

the variability of turbulence parameters in the marine boundary layer to 

establish a feeling for typical "background’* behavior. Nominal variations 
2 

of C^, and e during a 24 hour period (Fig. 3a and Fig. 3b) in the open ocean 

3 

are relatively smooth. Since e is proportional to u , we expect it to be 

2 

more variable than . A typical daily variation of wind velocity at 

2 

z = 10 meters ( u -j_q) is shown in Fig. 3c. Of course, is quite sensitive 

to changes in atmospheric temperature gradients as we can see in a com- 
parison of measurements in Monterey Bay under neutral conditions (7/76) 

2 

and unstable conditions (12/76) in Fig. 4. The diurnal variation of 

overland is dominated by solar heating of the surface during the day and 

2 

radiative cooling at night. This leads to a maximum of in the afternoon 

and minimum around sunrise and sunset. Since the ocean has both high heat 
capacity and excellent thermal conductivity (due to turbulent mixing), 
the surface temperature is subject to only a few tenths of degrees modu- 
lation by solar heating cycles. In the slightly unstable (R^ * - 0.1) 

boundary layer we typically encountered in the open ocean, we found that 

2 

radiation cooling of the atmosphere produced a slight maximum of during 

the night (Fig. 5a). In coastal areas, influenced by land-sea breeze cycles, 
we found diurnal variations more characteristic (but less extreme) of those 
found over land (Fig. 5b). 

B. Fog Data 

This paper will concentrate on the first four fog events encountered 
during CEVCOM 76: 

Fog //I Sept. 27 1100 - 0800 (Sept. 28) 
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Fog 


n 


0c t . 5 


0700 - 


1100 


Fog 


It 3 


Oct. 8 


1000 - 


1800 


Fog 


#4 


Oct. 9 


0800 - 


1200 



Fog event data are provided in tables in Appendix B. The location of the 

ACANIA during these events can be established from Fig. 6a (Fog //l) and 

Fig. 6b (Fog It's 2,3, and 4). Fog # 1 s 1,3, and 4 were open ocean fogs, 

all located west of Point Conception; Fog // 2 was a coastal fog located 

just south of Palos Verdes. Sea surface temperature contours (Fig. 7) for 

the first week of October were measured by Ralph Markson in the Airborne 

Research Associate's aircraft using a Barnes PRT-5 infrared radiometer 

(Backes, 1977). In order to simplify the presentation of four levels of 

turbulence data, the values from each level have been averaged together 

after translation to a 10 meter equivalent value using eqs. 10 and 12. 

The 10 meter equivalent averages of the temperature structure function 
2 

(<C T > ) for the four fog events are given in Fig. 8. Note the consider- 

10 2 

able variability of compared to Fig. 3a and Fig. 4. Figure 9 is a 

similar presentation for <£>^q. In some cases values of £ are not avail- 
able during fog since the impact of water droplets causes a noise spike 
by cooling the film. In heavy fog, this effect totally masks the real 
signal. We were able to obtain £ values during Fog // f s 1 and 4 by using 
hot wires which, due to their much smaller size (4.5 y diameter), are 
less affected. The heat fluxes for the fog events (Fig.s 10a and 10b) 
were calculated from the 10 meter averages using eqs. 15 and 18. Rep- 
resentative radiosonde dew point and temperature profiles are given in 
Fig. 11. 

The correlation of changes in turbulence parameters and fog events 
can be highlighted by preparing "ensemble 7 averages of the fog encounters. 
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Each fog event was divided into four one hour sections: one before and 

one after the fog is first encountered, and one before and one after the 

fog clears (Fog // 2 was excluded from this analysis). If a given fog event 

was less than two hours duration, the second and third sections were shortened 

2 

to 1/2 the fog duration. Note that both and e rise to a maximum about 

30 minutes before fog is encountered (Figs. 12a and 12b). One can also 

2 

examine the more fundamental behavior of in terms of fog and non-fog 

comparisons of height dependence (eq. 10) and the dimensionless temperature 
structure parameter (DTSP), f 2 (R i ), from eq . 11. Although the height de- 
pendence exhibits no obvious disimilarity (Fig. 13), the DTSP curve (Fig. 14) 

2 

implies lower values of for a given (3T/3z) during fog events. This 

is consistent with the lower fluxes obtained during fog by the turbulence 
method as compared to the gradient method. It is quite conceivable that 
this is an instrumental effect, caused by the accumulation of liquid water 
on the microthermal sensors. Such an accumulation, if large enough, could 
reduce the frequency response of the sensors sufficiently to attenuate the 
measured fluctuations. Based on laboratory measurements of frequency re- 
sponse in salt coated wires (Fairall & Schacher, 1977) one would expect 

this effect to be small until the water film was on the order of 30 microns 

2 

thick. The rapid increase of before dissipation of fog and the rapid 

decrease after dissipation are in the opposite direction of an expected 
water film effect. Despite the scarcity of droplet impacts observed on 
the 130 micron diameter hot films (implying the 2.5 micron diameter micro- 
thermal sensors remained very dry) during Fog # 2, a considerable drop in 
C was measured. 

j. 

Fog event data from other investigators can be found in Backes (1977) 
and Mack, et al . (1977). 
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V. Conclusions 



Despite the efforts of a number of scientists during a month of ship- 
board data gathering, unequivocal statements about the nature of marine 
fog are still difficult to come by. Although it might be wise to resist 
the temptation to embark into arm waving arguments about the reasons for 
fog formation, it is interesting to note that the three open ocean fogs 
described ‘were found in the vicinity of the cold water intrusion at Point 
Conception (upper left hand corner of Fig. 6) . To describe these fogs 
as being formed by "warm air blowing over cold water" is probably much 
too optimistic, since all were characterized by unstable surface layers. 
Indeed, the vertical heat fluxes often dramatically increased during the 
fog. We can make some statements about the nature of turbulence in fog, 
keeping in mind that they are subject to questions about the influence 
of water droplets on sensors. The temperature structure function (and 
therefore the turbulent heat transport) was subject to considerably more 

variations in the atmospheric conditions conducive to fog formation than 

2 

in the typical non-fog background. An ensemble analysis of and £ 

showed increased turbulence levels just before and just after the fog was 

encountered. This suggests that the fog is often superimposed on some 

large scale cellular structure in the boundary layer (on the order of 10 

km in the down wind direction) . When analyzed in dimensionless form 

2 

(normalized by temperature gradient), is lower in fog than out of fog. 

It is possible the fog droplets are retarding the temperature fluctuations 
by evaporation-condensation response. 

In general, we have made considerable progress towards the specific 
goals enumerated in the introduction. Notably in the provision of a 
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body of data for modelers to work with. However, there is considerable 
room for improvement (particularly in the water vapor flux capacity) in 
future work . 
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Appendix A 



The farms of the mean profile functions are taken from Businger, et al (1971). 



<P u (z/L) = (1 - 15 z/L) 1/4 z/L < o 

(j> u (z/L) = (1 + 4.7 z/L) z/L > o 

4> t (z/L) = .74 (1-9 z/L) _1/2 z/L < o 

4> t (z/L) = (.74 + 4.7 z/L) z/L > o 



The dimensionless velocity dissipation function is from Wyngaard and Cote 

(1971). 

E u (z/L) = [1 + .5|z/L| 2/3 ] 3/2 z/L < o 

E u (z/L) = [1 + 2.5(z/L) 3/5 ] 3/2 z/L > o 

The dimensionless temperature structure functions is from Wyngaard, et al 

(1971). 

f 1 (z/L) = 4.9 (1 - 7 z/L)' 2/3 z/L < o 

f-^z/L) = 4.9 (1 + 2.75 z/L) z/L > o 
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